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Abstract 

This  is  the  final  report  for  a  1-year  project  funded  by  the  European  Office  of  Research  and 
Development  (BOARD).  In  this  work  GaN/AlGaN  strain-balanced  heterostructures  have  been 
designed  for  near-infrared  absorption  through  intersubband  transitions.  By  making  use  of  existing 
collaborations  with  workers  at  the  University  of  Sheffield  and  UMIST  the  semiconductor  samples 
have  been  grown,  characterised  and  fabricated  into  devices.  The  experimental  results  have  been 
used  to  validate  the  design  procedure  and  the  prototype  devices  resulting  from  these  collaborations 
are  the  first  quantum  well  infrared  photodetectors  based  in  the  GaN  material  system  to  be  reported. 
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1  Summary  of  scientific  achievements 

1 .  Developed  a  model  for  calculating  the  conduction  band  potential  profi  le  of  a  GaN/Gai_ji:Al;tN 
heterostructure  including  the  spontaneous  polarisation  and  piezoelectric  h  elds[l] 

2.  Applied  the  model  to  the  design  of  quantum  wells  for  use  in  intersubband  devices  at  near- 
infrared  wavelengths  [  1  ] 

3.  Designed  quantum  well  systems  for  near-infrared  absorption  at  2,  3,  4  and  5  pm  for  use  in 
model  validation  experiments  [2] 

4.  Quantum  well  samples  were  grown  by  chemical  vapour  deposition  by  our  collaborator  Pe¬ 
ter  Parbrook  at  the  University  of  Sheffi  eld  and  characterised  with  Raman  spectroscopy  and 
photoluminescence  by  our  collaborator  Matthew  Halsall  at  UMIST[2] 

5.  Quantum  well  infrared  photodetectors  (QWIPs)  were  fabricated  at  UMIST  and  photocurrent 
measurements  were  performed 

6.  A  model  of  the  dark  (noise)  current  based  on  full  self-consistent  solution  of  the  subband  rate 
equations  was  developed  for  QWIPs.  This  was  validated  by  comparison  with  experiment[3] 


2  Scientific  Conclusions 

A  model  for  the  electron  energy  levels  in  strain-balanced  GaN/AlGaN  heterostructures  has  been  de¬ 
veloped  and  validated  by  comparison  with  experiment.  The  model  has  been  used  to  design  quantum 
well  structures  for  photoabsorption  in  the  near-infrared  spectral  range.  With  the  aid  of  collaborators 
these  semiconductor  heterostructures  have  been  grown  and  fabricated  into  the  h  rst  GaN  QWIPs.  In 
addition,  the  h  rst  fully  dynamical  model  of  the  dark  and  light  currents  in  QWIPs  has  been  developed 
and  applied  to  the  calculation  of  the  noise  current  and  responsivity.  The  results  of  these  calculations 
have  also  been  validated  by  experiment. 

3  List  of  publications  (attached) 

1.  V.  D.  Jovanonic,  Z.  Ikonic,  D.  Indjin,  P.  Harrison,  V.  Milanovic,  and  R.  A.  Soref,  Designing 
strain-balanced  GaN/AlGaN  quantum  well  structures:  Application  to  intersubband  devices  at 
1.3  and  1.55  pm,  J.  Appl.  Phys.,  93:3194  (2003) 

2.  M.  P.  Halsall,  B.  Sherliker,  P.  Harrison,  V.  D.  Jovanovic,  D.  Indjin,  Z.  Ikonic,  T.  Wang,  M.  A. 
Whitehead,  and  P.  J.  Parbrook,  Electronic  Raman  scattering  from  intersubband  transitions  in 
GaN/AlGaN  quantum  wells,  stat.  sol.  (c),  0:2662-2665  (2003) 

3.  V.  Jovanovic,  D.  Indjin,  Z.  Ikonic,  P.  Harrison  and  R.  A.  Soref,  Simulation  and  design  of 
GaN -based  near  infrared  quantum  well  infrared  photodetectors  ( QWIPs),  in  draft  form 
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List  of  conference  presentations  and  seminars 

1.  Condensed  Matter  and  Materials  Physics,  Belfast,  April  2003,  ‘Design  issues  for  GaN 
intersubband  devices’,  presented  by  V.  Jovanovic 

2.  7th  International  Conference  on  Intersuhhand  Transitions  in  Quantum  Wells,  Switzer¬ 
land,  September  2003,  ‘X  ~  3.5pm  intersubband  transitions  and  electronic  Raman  scattering 
study  in  strain  balanced  GaN/AlGaN quantum  wells’,  presented  by  V.  Jovanovic 

3.  U.K.  Nitrides  Consortium  Meeting,  Cardiff,  June  2003,  ‘Recent  progress  in  GaN/AlGaN 
intersubband  devices  design  ’ 


5  Suggestions  for  future  projects 

This  project  has  been  a  feasibility  study  to  see  if  GaN-based  QWIPs  with  suffi  ciently  large  conduc¬ 
tion  band  offset  could  be  designed  for  near-infrared  absorption.  The  scientifi  c  output  of  the  project 
has  demonstrated  that  this  is  the  case  and  this  has  been  further  substantiated  with  the  growth  and 
characterisation  of  several  samples.  This  work  is  now  at  a  crossroads  and  as  this  report  is  written  no 
funds  have  been  secured  with  which  to  continue  this  unique  and  leading  work. 

The  opportunity  is  therefore  taken  in  this  fi  nal  report  to  suggest  potential  projects  and  the  re¬ 
quired  level  of  funding  in  order  to  meet  certain  targets: 

1.  Optimised  uncooled  near-infrared  single-pixel  detectors  operating  across  the  1-5  pm 
wavelengths  In  this  project  the  fi  rst  GaN-based  QWIPs  have  been  designed,  fabricated  and 
characterised.  The  short  duration  of  the  project  (1  year)  and  the  limited  budget  have  only 
allowed  for  one  iteration  in  the  design  and  fabrication  cycle  and  further  improvements  in  the 
performance  of  the  devices  would  require  further  time  and  funding.  A  three  year  programme 
with  3  PhD  students  would  be  an  effective  way  to  continue  in  this  line  of  research  and  it  is 
likely  that  the  objective  of  uncooled  detection  of  near-infrared  radiation  could  be  achieved 
with  a  single  pixel  device. 

Funding  required:  Around  $100k/year  for  3  years 

2.  Focal  plane  arrays  (FPAs)  for  near-infrared  imaging  and  spectroscopy  This  more  ambi¬ 
tious  aim  would  be  attempted  as  a  follow  on  from  the  above.  This  would  require  the  optimi¬ 
sation  of  the  single  pixel  detectors  and  then  the  development  of  the  fabrication  technologies 
to  achieve  large  scale  (10,000-100,000)  arrays  of  pixels  for  connection  to  commercially  avail¬ 
able  readout  circuits  and  cameras.  This  has  been  achieved  in  a  selection  of  other  materials 
including  GaAs-  and  HgCdTe-based  devices.  The  resources  required  for  this  are  likely  to  be 
2  Postdoctoral  level  researchers  for  3  years.  Funding  required:  Around  $150k/year  for  3 
years 
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A  criterion  for  strain  balancing  of  wurtzite  group-III  nitride-based  multilayer  heterostructures  is 
presented.  Single  and  double  strain-balanced  GaN/AIGaN  quantum  well  structures  are  considered 
with  regard  to  their  potential  application  in  optoelectronic  devices  working  at  communication 
wavelengths.  The  results  for  realizable,  strain-balanced  structures  are  presented  in  the  form  of 
design  diagrams  that  give  both  the  intersubband  transition  energies  and  the  dipole  matrix  elements 
in  terms  of  the  structural  parameters.  The  optimal  parameters  for  structures  operating  at  \  ~  1 .3  and 
1.55  jjm  were  extracted  and  a  basic  proposal  is  given  for  a  three  level  intersubband  laser  system 
emitting  at  1.55/U.m  and  depopulating  via  resonant  longitudinal  optical  (LO)  phonons  (^(Ulo”“90 
meV).  ©  2003  American  Institute  of  Physics.  [DOT  10.1063/1.1556177] 


I.  INTRODUCTION 

In  recent  years  intersubband  transitions  (ISBTs)  in 
AlGaN-based  heterostructures  have  emerged  as  a  promising 
route  for  access  to  the  near-infrared  spectral  range,  while 
also  showing  ultrafast  carrier  dynamics.  A  few  groups  have 
reported  intersubband  absorption  for  wavelengths  around  1 .3 
and  1.55  /rm,*“^  that  are  important  for  optoelectronic  tele¬ 
communications,  as  well  as  ultrashort  relaxation  times 
~fs  in  quantum  well  structures  and  superlattices.*  These 
experimental  successes  suggest  possible  application  in  quan¬ 
tum  well  infrared  photodetectors  (QWIPs),  optical  switches 
or  even  intersubband  lasers.  Improvements  in  the  molecular 
beam  epitaxy  growth  technique  that  allow  control  of  layer 
widths  of  the  order  of  a  few  monolayers  make  feasible  the 
design  and  growth  of  complex  structures.  However,  growth 
of  uniform,  dislocation-free  multilayer,  optically  thick  stacks 
is  still  a  considerable  challenge."  *^  To  achieve  pseudomor- 
phic  growth,  there  is  the  requirement  of  strain  balance  upon 
the  whole  structure.  Although  multilayer  structures  generally 
show  much  better  stability  against  relaxation  than  do  homo¬ 
geneous  strained  layers,  one  should  still  aim  at  achieving  full 
strain  balance  when  designing  such  a  structure  in  order  to 
have  reliable,  reproducible  growth. 

Several  approaches  to  the  design  of  strain-balanced 
structures  have  been  considered,**  but  they  were  usually  re¬ 
stricted  to  two-layer  structures  (e.g.,  single  quantum  wells) 
based  upon  semiconductor  materials  with  cubic  symmetry.  In 
this  article  we  extend  the  “zero-strain”  approach**  to 
multilayer  structures  based  upon  wurtzite  (e.g.,  GaN/AlN) 
semiconductors  grown  in  the  conventional  [0001]  direction. 


“'Electronic  mail:  eenvj@leeds.ac.uk 


and  then  we  present  design  diagrams  oriented  their  possible 
application  in  optoelectronics. 


II.  THEORETICAL  CONSIDERATIONS 

The  zero-strain  condition**  is  the  only  theoretical  ap¬ 
proach  that  accounts  for  different  elastic  constants  of  layers, 
in  contrast  to  the  more  common  average  lattice  or  thickness 
weighted  methods.*"*  It  is  based  on  minimizing  the  average 
elastic  energy  for  the  whole  multilayer  structure,  i.e.. 


(1) 


where  is  the  width  and  /j.  strain  energy  density  of  the  kth 
layer. 


Within  the  framework  of  classical  elastic  theory 
elastic  energy  density  of  the  kth  layer  is 
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where  c]*'  is  the  elastic  tensor  element  and  the  strain 
tensor  element  for  a  hexagonal  crystallographic  system  in 
matrix  notation**  of  the  kth  layer. 

Since  we  consider  a  wurtzite  based  multilayer  system 
grown  along  the  [0001]  axis,  the  corresponding  stress  tensor 
element  (0-3)  should  equal  zero.  By  using  formulas  given  in 
Refs.  15  and  16  the  biaxial  strain  tensor  has  only  diagonal 
elements  and  is  given  by 
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where  ei  =  (a  — ao)/ao  the  in-plane  strain,  a  is  the  buffer 
in-plane  lattice  constant  (along  the  wurtzite  a  axis)  and  Og  is 
the  lattice  constant  of  a  layer  in  strain-free  conditions.  From 
Eqs.  (2)  and  (3)  it  follows  that  the  strain  energy  density  in  a 
layer  is 

(4) 

where  A  =  Cii  +  Ci2~2c[3/c33 .  Equation  (1)  may  then  be 
written  as 

J  -srn  j  ‘  / 

^k=l^k 

The  pseudomorphic  condition  relates  the  strain  in  adjacent 
layers  (e®  and  e**  and  is  given  as 


and  the  average  in-plane  stress  that  follows  is 
-  df  2 
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The  zero-stress  (or  strain-balance)  condition  means  that  the 
in-plane  stress  is  zero,  and  using  =  ui! a^.  (k 

=  2,3,  .  .  . ,)  Eq.  (7)  then  gives 


A,e[%+A,e[^\"-^+AA%-  +  ,--,  =  0,  (8) 

G2  ^3 

which  delivers  the  buffer  lattice  constant  necessary  for  the 
strain-balanced  multilayer  structure  as 

^k=iAklk/ak 

ao  = - (9) 

n=,Aklklal 

and,  therefore,  in  case  of  GaN/AlGaN  structures,  the  Al  con¬ 
tent  in  the  buffer  follows  straightforwardly. 

If  a  multilayer  structure  has  periodicity,  the  above  ex¬ 
pression  should  be  understood  to  apply  to  n  layers  that  con¬ 
stitute  a  single  period,  and  the  whole  structure  will  then 
clearly  be  strain  balanced. 

The  strain  level  in  individual  layers  directly  influences 
the  conduction  band  profile,  due  to  the  large  strain- 
dependent,  piezoelectric  polarization  induced  electric  field. 
This  means  that  strain-balanced  device  design  and  optimiza¬ 
tion  are  more  complicated  than  would  be  the  case  with  no 
polarization.  This  is  because  any  single  structural  parameter 
cannot  be  varied  independently  (in  order  to  change  the  sub¬ 
band  structure),  so  the  strain-balance  condition  has  to  be  met 
as  well.  In  the  remaining  part  of  this  article,  we  present  suit¬ 
able  design  diagrams  for  strain-balanced  single  and  double 
GaN/AlGaN  QW  structures;  see  Eig.  1 .  Single  QW  structures 
were  explored  with  respect  to  the  ISBT  energies  and  dipole 
matrix  elements  corresponding  to  the  1— >2  transition  with  a 
view  towards  their  possible  application  in  devices  like 
QWIPs  or  saturable  switches.  Double  QWs  were  explored 


FIG.  1.  Conduction  band  profile  and  wave  function  moduli  of  single  and 
coupled  double  GaN/AlGaN  QWs,  and  the  notation  of  structural  parameters. 


with  respect  to  transition  energies  between  |l)  and  |2),  and 
1 2)  and  |3),  with  a  view  towards  their  potential  application 
in  three-level  intersubband  lasers.** 

In  the  calculations  we  have  employed  the  effective  mass 
model  and  solved  the  envelope  function  Schrbdinger  equa¬ 
tion  with  bulk  nonparabolicity  included  via  the  energy  de¬ 
pendent  effective  mass  (two  band  Kane  model). *^  Nonlinear 
piezoelectric  and  spontaneous  polarizations  were  also 
included^**  and  the  fields  induced  were  calculated  from  the 
simple  expression^*  (all  constraints  of  that  approach  are  sat¬ 
isfied  here^^) 


p  _^k{Pk  P j)hl^k 


(10) 


where  P ,  =  is  the  total  polarization  given  as  a  sum 

of  piezoelectric  and  spontaneous  (pyroelectric)  polarization 
within  a  layer  of  width  1, .  In  calculating  the  alloy  band  gap 
we  have  included  the  effect  of  bowing.^* 

The  Schrbdinger  equation  used  within  this  model  is 
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+  Uiz)'i'  =  E'if, 
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where  the  position-  and  energy-dependent  effective  mass  is 
m(z,E)  =  m(z){l+[E-U(z)yEg(z)},  where  m(z),  Eg(z) 
and  U(z)  denote,  respectively,  the  parabolic  effective  mass, 
band  gap  and  conduction  band  edge  (which  includes  the 
band  edge  discontinuities  and  the  polarization-induced  po¬ 
tential).  All  three  vary  along  coordinate  z  because  of  modu¬ 
lation  of  the  material  composition.  This  equation  may  be 
solved  by  a  standard  shooting  method,  starting  with  a  well- 
behaved  wave  function  on  one  side  and  scanning  for  energies 
where  it  is  well  behaved  on  the  other  side.  However,  this 
method  loses  stability  for  wider  structures,  typically  with 
those  five  wells  or  so.  We  have  therefore  employed  a  finite- 
difference  method,  in  which  Eq.  (11)  is  written  in  discretized 
form,  in  points  along  the  z  axis  like,  e.g.,  in  Ref.  24,  and 
reduces  to  a  matrix  eigenvalue  problem,  where  the  Hamil¬ 
tonian  matrix  is  tridiagonal.  However,  the  energy  E  appears 
not  only  in  the  main  diagonal,  but  in  all  the  terms  (via  the 
effective  mass),  which  leads  to  a  nonlinear  eigenvalue  prob¬ 
lem. 


[H{E)]i/k=Eif,, 


(12) 
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which  cannot  be  solved  by  conventional  diagonalization  rou¬ 
tines  that  handle  linear  eigenvalue  problems  [//]i/r=£i/' 
(which  would  be  the  case  if  the  Hamiltonian,  i.e.,  the  effec¬ 
tive  mass  in  it,  is  not  energy  dependent).  In  the  field  of  elec¬ 
tronic  structure  calculations  nonlinear  eigenvalue  problems 
also  appear,  e.g.,  in  augmented  plane  wave  methods,  but  the 
matrix  there  is  not  tridiagonal,  and  the  computational  de¬ 
mands  are  considerably  larger  than  for  linear  eigenproblems, 
so  sophisticated  methods  have  been  devised  to  speed  up  the 
process. Here,  however,  we  take  advantage  of  the  tridiago¬ 
nal  form  of  the  finite-difference  representation  Eq.  (11)  for 
fast  solving  of  Eq.  (12).  We  first  scan  for  energies  £,  where 

det([//(£,.)] -£,[/])  =  0  (13) 

is  satisfied  [the  £,  are  obviously  the  eigenvalues  of  Eq.  (12) 
and  [/]  is  the  unity  matrix].  Calculation  of  the  deter¬ 

minant  of  the  tridiagonal  matrix^®  is  fast  and  scales  linearly 
with  its  size  (just  as  does  the  shooting  method),  hence  it  is 
easy  to  achieve  high  precision.  In  this  particular  problem  it 
suffices  to  find  just  the  sign,  not  the  full  value  of  the  deter¬ 
minant,  and  supply  it  to  the  bisection  routine,  a  feature  which 
almost  cuts  into  half  the  computational  time.  Normally,  one 
is  interested  in  only  a  few  lowest,  rather  than  all,  eigen¬ 
values  of  Eq.  (12). 

To  find  the  eigenvectors,  i.e.,  the  wave  functions,  Eq. 
(12)  can  be  written,  for  any  particular  £,  in  the  form  of  a 
fictitious  linear  eigenproblem, 

([//(£,.)]-£,.[/])  (A,- =  X*A,--  (14) 

One  of  the  eigenvalues  X.  will  be  (almost,  because  of 
roundoff  errors)  equal  to  zero,  and  the  eigenvector  corre¬ 
sponding  to  it  is  actually  the  eigenvector  of  Eq.  (12)  corre¬ 
sponding  to  eigenenergy  £,  [other  eigenvalues  and  eigenvec¬ 
tors  of  Eq.  (14)  are  physically  meaningless].  Eor  this  task  we 
use  standard  diagonalization  routines  (EISPACK)  that  allow 
evaluation  of  a  single,  selected  eigenvector — that  with  the 
already  known,  zero  eigenvalue.  Alternatively,  upon  finding 
eigenvalues  £, ,  one  may  also  calculate  eigenvectors  by  solv¬ 
ing  the  tridiagonal  system  of  linear  equations  corresponding 
to  Eq.  (12)  using  appropriate  fast  routines,^®  but  this  ap¬ 
proach  was  found  to  offer  somewhat  less  stability. 

The  method  was  found  to  be  about  as  equally  fast  and 
accurate  as  the  shooting  method  (in  cases  where  both  work 
well),  while  retaining  full  stability  where  the  shooting 
method  fails. 

The  calculated  results  were  used  to  make  contour  plots, 
from  which  the  sets  of  structural  parameters  that  deliver  the 
transition  energy  or  matrix  element  of  interest  may  be  easily 
read,  thus  allowing  some  optimal  choice  to  be  made. 


III.  NUMERICAL  RESULTS  AND  DISCUSSION 

Eigure  2  shows  contour  plots  of  both  the  ISBT  energies 
and  dipole  matrix  elements  squared,  corresponding  to  the  1 
^2  transition  in  a  single  QW  structure  as  functions  of  the 
relevant  structural  parameters  [for  the  notation  see  Eig.  1  (a)], 
i.e.,  the  well  width  (l„)  and  the  A1  content  in  the  barrier 
layer  (xh’^Ui,/ AE^ ,  where  AE^  is  the  conduction  band  off- 
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FIG.  2.  Contour  plot  of  the  transition  energies  (solid  line)  and  dipole  matrix 
element  squared  (dashed  line)  between  levels  |l)  and  |2)  in  single  AlGaN/ 
GaN  QWs  obtained  for  a  range  of  structural  parameters. 

set  set  to  2  eV).^^^*  ^  fixed  value  of  the  barrier  width  was 
chosen  as  Z;,  =  60  A  since  it  does  not  influence  the  results 
much. 

In  order  to  give  specific  examples,  we  consider  the  cases 
of  1.3  (  —  950  meV)  and  1.55/U,m  (  —  800  meV)  transitions, 
which  are  important  in  applications.  The  structural  param¬ 
eters  that  deliver  the  1.55/U,m  intersubband  absorption  can  be 
read  from  the  solid  line  labeled  800.0  (in  units  of  meV)  in 
Eig.  2.  Any  choice  along  this  line  gives  the  same  wavelength, 
but  the  dipole  matrix  element  di2  varies.  There  are  two  sets 
giving,  e.g.,  di2=  16  A  ;  these  are  at  the  intersections  of  the 
dashed  line  labeled  16  with  the  800.0  solid  line.  Usually  one 
aims  to  have  the  largest  possible  di2,  and  these  (optimal) 
structural  parameters  occur  where  a  dashed  line  just  touches 
the  800.0  solid  line;  in  this  case  it  is  the  dashed  line  labeled 
20,  i.e.,  the  maximal  value  of  c/[2=20  A  and  it  is  obtained 
with  l„=  18  A  and  Xi,  =  Q.S  (this  point  is  marked  by  a  circle 
and  an  arrow  in  Eig.  2).  Similarly,  one  can  see  from  Eig.  2 
that  the  absorption  at  1.3yU,m  can  be  achieved  only  in  a  rather 
narrow  range  of  parameters:  17A  and  X|^>0.85,  and  the 

largest  possible  value  of  t/i2=  16  A  is  obtained  for  l„=  17 
A  and  x^  =  0.94  (also  marked  in  Eig.  2).  Along  with  the  op¬ 
timal  parameters,  Eig.  2  shows  the  range  of  parameters  that 
can  provide  any  (suboptimal)  structure  for  a  given  wave¬ 
length,  and  also  indicates  where  the  sensitivity  to  a  param¬ 
eter  becomes  only  marginal.  These  may  be  important  when 
the  optimal  design  cannot  be  fabricated,  for  whatever  reason. 

The  contour  plots  shown  in  Eigs.  3  and  4  show  ISBT 
energies  corresponding  to  the  1-^2  and  2^3  transitions  in 
double  QW  structures  as  a  function  of  their  structural  param¬ 
eters  [for  the  notation  see  Eig.  1(b)],  i.e.,  the  first  well  (l„i) 
and  middle  barrier  width  (li,),  for  two  fixed  values  of  the 
second  well  width  (l„2~^  8  A).  In  this  calculation  the 

outer  barriers  were  chosen  to  be  of  AIN,  in  order  to  accom¬ 
modate  the  three  states,  and  their  widths  were  also  set  to  60 
A.  An  interesting  case,  in  view  of  the  possible  application  of 
double  QW  systems  as  active  regions  of  intersubband  lasers, 
occurs  when  the  energy  of  the  1^2  transition  is  close  to  the 
longitudinal  optical  (LO)  phonon  energy  (^&)lo^90  meV  in 
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FIG.  3.  Contour  plot  of  transition  energies  between  |l)  and  |2)  (dashed 
line)  and  |2)  and  |3)  (solid  line)  states  in  double  AlN/GaN  QWs  obtained 
for  a  range  of  structural  parameters.  The  second  well  width  is  set  to  l„2 
=  6  A. 


GaN).  This  provides  ultrafast  depopulation  of  the  second 
subband  and  the  possibility  of  population  inversion  between 
levels  1 3)  and  |2).  An  interesting  choice  is  if  the  spacing 
between  the  second  and  third  subband  is  set  to  a  wavelength 
of  1.55/U.m,  which  is  important  in  fiber-based  communica¬ 
tions.  For  l„2~^  A  this  situation  appears  for  the  first  well 
width  of  =  19  A  and  the  middle  barrier  width  of  17  A.  A 
wider  second  well  A)  requires  a  much  wider  first 

well,  l„i  =  29  A,  while  the  middle  barrier  width  is  1/,=  16  A. 
That  is  a  direct  consequence  of  the  presence  of  polarization 
induced  fields  in  the  layers  and  the  therefore  inherent  tilt  of 
the  structure. 

For  any  set  of  parameters  from  these  diagrams,  the 
strain-balance  requirement  can  be  straightforwardly  deter¬ 
mined  from  Eq.  (9). 
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FIG.  4.  Same  as  in  Fig.  3,  but  for  lw2~^ 


IV.  CONCLUSION 

A  strain-balance  condition  for  wurtzite  GaN-based 
multilayer  structures  was  been  derived  and  useful,  “design- 
friendly”  contour  plots  for  single  and  double  QW  structures 
complying  with  the  strain-balance  requirement  are  given. 
From  these  plots  optimal  parameters  were  extracted  for  ISBT 
absorption-based  devices,  as  well  as  for  a  candidate  structure 
for  a  three  level  intersubband  laser. 
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We  present  a  Raman  scattering  study  of  two  GaN/Al^Ga,_^N  multiple  quantum  well  (MQW)  structures 
with  X  =  0.3  and  x  =  0.4  and  well  widths  of  6  nm  and  4  nm  respectively,  they  were  nominally  undoped  but 
are  expected  to  contain  a  low-density,  n-type,  earrier  population  due  to  the  residual  donors  in  the  barriers. 
Polarization  dependent  Raman  scattering  was  performed  at  room  temperature,  strong  scattering  due  to  in¬ 
tersubband  transitions  (ISBT’s)  in  the  GaN  quantum  wells  was  observed  from  both  the  ei-e2  and  ej-es 
transitions  for  both  samples.  The  first  sample  has  an  ei-e2  transition  at  300  meV  and  ei-e3  at  480  meV 
whilst  the  second  sample  has  an  ei-e2  transition  at  360  meV  and  ej-es  at  580  meV.  The  fall  width  half 
maximum  of  the  Raman  peaks  is  lower  than  that  reported  in  infra-red  absorption  for  similar,  heavily  do¬ 
ped,  MQW  structures. 
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1  Introduction  Electronic  devices  that  utilize  intersubband  transitions  (ISBT’s)  are  of  increasing  tech¬ 
nological  importance  for  applications  in  the  Far  and  Mid  infrared  regions  of  the  spectrum.  The  emis¬ 
sion/absorption  wavelength  of  an  ISBT  is  only  limited  by  the  size  of  the  conduction  band  offsets  of  the 
material  system  used.  Moreover,  The  devices  are  unipolar  and  therefore  only  the  control  of  n-type  dop¬ 
ing  is  required.  Since  the  first  demonstration  of  an  intersubband  laser  (ISBL)  in  the  InGaAs/GaAs  system 
a  number  of  other  111-V  systems  have  been  used  to  push  the  wavelength  down  to  around  3  pm  [1-3].  In 
this  respect  the  AlGaN/GaN  system  is  the  center  of  much  current  interest  as  the  large  offsets  in  this  sys¬ 
tem  allow  the  production  of  ISBT’s  corresponding  to  wavelengths  as  short  as  1.35  pm  [4].  An  intersub¬ 
band  laser  based  on  such  technology  would  in,  addition  to  having  wide  wavelength  tunability  and  high 
optical  power,  benefit  from  the  fast  carrier  dynamics  allowing  terahertz  modulation  rates  [5].  Although 
ISBT’s  have  been  reported  in  this  region  of  the  spectrum,  the  production  of  such  devices  has  not  been 
demonstrated.  This  is  in  large  part  due  to  the  difficulties  in  operating  with  a  system  with  such  large  lat¬ 
tice  mismatches  and  with  the  thin  layers  needed  to  achieve  short  wavelength  ISBT’s.  For  progress  to  be 
made  in  the  system  prospective  device  structures  need  to  have  their  electronic  transitions  characterized  as 
completely  as  possible. 

The  usual  technique  used  to  study  ISBT’s  is  infrared  absorption.  Experimentally  the  technique  is  rela¬ 
tively  straightforward  and  is  generally  either  done  by  turning  the  sample  to  the  Brewster  angle  and  meas¬ 
uring  the  p-polarised  transmission  or  by  waveguiding  the  light  into  the  sample.  The  former  technique 
works  well  in  the  case  of  highly  doped  sample  with  large  numbers  of  quantum  well  repeats,  so  that  sig- 
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nificant  absorption  can  be  measured  [6].  The  use  of  waveguiding  is  more  problematic  experimentally  but 
is  capable  of  measuring  small  absorption  changes.  Indeed,  absorption  due  to  61-62  IS  transition  in  a 
nominally  undoped  AlGaN/GaN  high  eleetron  mobility  structure  (single  quantum  well)  were  recently 
reported  in  this  manner  [7]  as  well  as  absorption  due  to  the  61-63  transition  in  MQW  structure  [8]. 

Here,  we  report  an  electronic  Raman  scattering  (ERS)  study  of  two,  nominally  undoped,  GaN/AlGaN 
multiple  quantum  well  (MQW)  samples  designed  to  have  ISBT’s  corresponding  to  wavelengths  in  the 
vicinity  of  3  and  4  pm  and  with  strain  balancing  condition  enforced  by  employing  zero-strain  approach 
[9].  The  nominally  undoped  samples  are  expected  to  contain  a  low  background  n-type  carrier  concentra¬ 
tion;  our  own  IR  absorption  measurements  were  unable  to  observe  the  ISBT’s  due  to  the  low  carrier 
concentrations.  The  use  of  the  ERS  technique,  however,  allows  the  observation  of  both  the  61-62  and 
61-63  transitions  of  the  structures  without  any  further  sample  preparation.  Moreover,  Raman  selection 
rules  allow  the  observation  of  even  parity  transitions  which  are  nominally  forbidden  in  absorption  meas¬ 
urements.  The  observed  transition  energies  are  modelled  using  an  effective  mass  model  combining  the 
effect  of  intrinsic  electric  fields  with  a  self-consistent  solution  to  Poisson’s  equation  for  an  assumed,  low, 
carrier  concentration. 

2  Experimental  procedure 

2.1  Sample  growth  All  investigated  samples  were  grown  on  (0001)  sapphire  substrates  by  low  pres¬ 
sure  metalorganic  vapour  phase  epitaxy.  Sapphire  substrates  were  initially  treated  in  H2  ambient  at 
1 140  °C,  followed  by  the  growth  of  a  25  nm  thick  low-temperature  (550  °C)  GaN  nucleation  layer  and  a 
1.5  pm  thick  layer  of  nominally  undoped  GaN  grown  at  1050  °C.  Prior  to  the  growth  of  the  AlGaN/GaN 
superlattice  structure  at  1050  °C,  an  AIN  interlayer  of  about  20  nm  was  grown  at  600  °C  in  order  to  pre¬ 
vent  the  crack  formation  due  to  the  large  lattice  and  thermal  expansion  coefficient  mismatch  between 
AlGaN  and  GaN.  The  MQW  structure  of  the  samples  grown  is  detailed  in  Table  1. 


Table  1  Sample  parameters. 


sample  no. 

no.  of  periods 

well  width 

barrier  width 

1 

10 

6  nm  GaN 

28  nm  Alo  sGao.TN 

2 

30 

4  nm  GaN 

6  nm  Al()  42Gao.58N 

2.2  Raman  technique  The  Raman  spectra  were  recorded  in  backscattering  geometry  from  the  cleaved 
sample  edge  using  a  Raman  microscope  system.  The  514  nm  line  of  an  argon-ion  laser  was  used  for 
excitation.  The  spectra  were  analysed  for  polarization  with  the  growth  direction  with  the  polarization 
axis  being  parallel  to  the  growth  direction  (x)  or  at  right  angles  to  it  (y).  The  spectra  were  recorded  at 
room  temperature. 

3  Results  and  discussion  Figure  1  shows  the  Raman  spectra  obtained  from  Sample  1,  the  upper  spec¬ 
trum  was  recorded  under  polarized  conditions  (using  standard  notation)  [z(y,  y)  z]  and  the  lower  spec¬ 
trum  under  depolarized  [z{x,y)z].  In  the  figure,  there  are  clearly  two  features  at  300  meV  and 
480  meV,  which  are  observed  under  depolarized  conditions,  as  is  appropriate  for  Raman  scattering  from 
excitations  of  single  carrier  character  [10,  11].  Figure  2  shows  an  equivalent  series  of  spectra  for  sample 
2.  The  peaks  at  500-600  meV  (Raman  shift)  in  both  spectra  are  actually  emission  lines  due  to  Chro¬ 
mium  impurities  in  the  sapphire  substrate.  In  the  figure  the  top  trace  is  again  recorded  in  the  depolarised 
geometery,  the  middle  trace  is  recorded  in  under  polarised  conditions,  the  bottom  trace  is  a  subtraction  of 
the  polarised  from  the  depolarised  scaled  so  as  to  remove  the  chromium  emission  from  the  spectrum. 

Again  in  the  depolarized  spectrum  we  see  a  feature,  which  this  time  occurs  at  360  meV.  In  the  sub¬ 
tracted  spectrum  there  is  also  a  feature  at  580  meV,  which  is  obscured  the  depolarised  spectrum  by  the 
chromium  emission.  The  observed  features  from  both  samples  lie  close  to  the  expected  values  for  the 
61-62  and  61-63  transitions  of  each  sample  respectively  and  we  ascribe  them  to  electronic  Raman  scatter- 
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Raman  shift  (meV) 

Fig.  1  Raman  spectrum  of  sample  1  recorded  in 
both  polarised  (upper  trace)  and  depolarised 
(lower  trace)  scattering  geometries. 


Fig.  2  Raman  spectrum  of  sample  2  recorded  in  depolar¬ 
ised  (upper  trace)  and  polarised  (middle  trace)  scattering 
geometries.  The  lower  trace  is  a  scaled  subtraction  of  the 
polarised  spectrum  from  the  depolarised. 


ing.  The  linewidths  of  these  features  are  remarkably  narrow,  for  example  in  sample  1  ei-e2  has  a  transi¬ 
tion  energy  of  300  meV  and  a  FWHM  of  20  meV.  In  reference  5  ISBT  absorption  was  reported  from 
ei-e2  at  306  meV  from  a  nominally  undoped  AlGaN/GaN  HEMT  structure  with  a  FWHM  of  86  meV, 
and  in  6  a  heavily  doped  MQW  structure  was  reported  to  show  ei-e2  absorption  at  300  meV  with  a 
FWHM  of  96  meV.  It  can  be  seen  in  the  figures  that  scattering  intensity  due  to  the  Ci-Cs  transition  is 
substantially  stronger  than  that  due  to  the  ei-e2  transition.  The  same  effect  has  been  observed  in  the 
modulation  doped  AlAs/GaAs  system  [11]  and  is  a  result  of  the  Raman  matrix  element  for  scattering  via 
spin  density  excitations  being  non-zero  only  for  even  parity  excitations  in  the  dipole  approximation.  The 
odd  parity  ei-e2  transition  is  only  allowed  via  high  order  quadropole  terms.  This  is  not  necessarily  true  in 
the  presence  of  strong  resonant  Raman  effects,  in  the  this  current  system  it  is  unclear  which  electronic 
states  could  be  resonant  with  2.42  eV  excitation  and  we  assume  that  the  scattering  is  non-resonant. 

The  observation  of  electronic  Raman  scattering  has  been  used  in  other  systems  to  obtain  accurate 
values  of  the  band  offsets  in  the  absence  of  uncertainties  due  to  exciton  binding  energies  etc.  The  ob¬ 
served  transition  energies  have  been  modelled  self-consistently  by  iterative  solution  of  the  Schrodinger 
and  Poisson  equations  with  effects  of  conduction  band  nonparabolicity  included  via  Kane  two  band 
model  of  an  energy  dependent  effective  mass  [12,  13].  The  large  built-in  electric  fields  (~MV/cm)  in¬ 
duced  by  nonlinear  piezoelectric  and  spontaneous  polarization  was  taken  into  account  and  calculated 
from  the  simple  expression  assuming  periodic  boundary  conditions  over  the  whole  structure  [14].  All 
ternary  alloy  constants  were  estimated  by  employing  Vegard’s  linear  law.  We  assume  a  conduction  band 
offset  Ub  and  that  Al  content  in  the  barrier  layer  is  given  as  x  =  UiJAE„  where  AE^  is  the  AlN/GaN  con¬ 
duction  band  offset,  set  to  be  ~2  eV.  In  calculating  the  alloy  band  gap  we  have  included  the  effect  of 
bowing  (ft  =  1  eV)  [15]. 

By  using  the  described  self-consistent  model  we  calculated  transition  energies  based  on  the  specified 
sample  structural  parameters  varied  within  the  growth  uncertainty.  The  calculated  values  of  Cj-e,  and 
ej-e,  transition  energy  are  308  meV  and  390  meV  for  sample  1  and  390  meV  and  490  meV  for  sample  2. 
It  can  be  concluded  that  estimated  values  for  ej-e,  are  in  very  good  agreement  with  measured  transition 
energies  of  an  average  5%  while  there  is  a  certain  discrepancy  (>10%)  observed  for  03-6,  transitions. 
This  is  probably  due  to  calculated  intrinsic  field  value  uncertainty  in  the  barrier  layer  which  for  lower  Al 
content  significantly  influence  the  energy  position  of  the  second  excited  state.  From  the  measured  data 
by  fitting  to  values  for  ej-e,  transition  energies,  conduction  band  offset  can  be  estimated  to  be  of  around 
1.95  eV  for  sample  2  and  around  1.8  eV  for  sample  1. 

In  conclusion,  we  have  demonstrated  the  application  of  Raman  scattering  to  the  study  of  nitride  semi¬ 
conductor  intersubband  transitions.  Two,  nominally  undoped,  structures  were  studied  and,  despite  the 
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low  residual  doping  level,  seattering  due  to  transitions  from  the  ground  state  to  the  first  two  excited  state 

electron  subbands  was  readily  observed  from  both  samples. 
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Abstract 

A  model  for  electron  transport  in  quantum  well  infrared  photodetectors  with  more  than  one  discrete 
state  is  presented,  based  on  a  scattering  rate  equation  analysis.  The  important  macroscopic  parameters 
like  current  density,  responsivity  and  capture  probability  can  be  estimated  directly  from  a  fi  rst  principals 
calculation.  The  applicability  of  the  model  was  tested  by  comparison  with  experimental  measurements 
from  a  GaN/AlGaN  device,  and  good  agreement  was  found.  The  model  in  general  can  be  applied  to  any 
other  material  system  or  QWIP  design. 
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INTRODUCTION 


The  recent  observations  of  near-infrared  intersubband  transitions  in  wide-band  gap  group-III 
nitride  semiconductor  quantum  wells  (QWs)  [1-6],  as  well  as  the  experimental  verifi  cation  of  the¬ 
oretically  predicted  femtosecond  carrier  dynamics  [7-9],  confirm  that  GaN-based  materials  offer 
a  route  for  the  design  of  complex  optoelectronic  devices  like  ultrafast  waveguide  switches,  opti¬ 
cal  modulators  and  even  quantum  cascade  lasers.  Moreover,  rapid  improvements  in  the  growth 
technology,  primarily  by  using  Molecular  Beam  Epitaxy,  make  possible  the  growth  of  multi¬ 
layer  structures  with  considerable  quality  of  the  epitaxial  layers,  delivering  controllable  widths 
of  the  order  of  a  monolayer  [10,  11].  This  enables  investigation  of  even  more  complicated  inter¬ 
subband  devices  [12],  like  quantum  well  infrared  photodetectors  (QWIPs),  which  would  both 
extend  the  present  range  of  materials  for  QWIPs  and  also  the  detection  wavelengths  towards 
the  optical  communications  range  around  1.55  /urn.  Previous  experience  in  the  design  of  novel 
quantum  devices  [13]  clearly  shows  that  systematic  theoretical  modelling  is  a  necessary  step  to¬ 
wards,  improvements  of  the  existing  structures  [14],  and  the  understanding  of  physical  processes 
within  [15,  16]. 

In  this  paper,  two  main  topics  are  addressed:  (i)  development  of  a  model  for  the  simulation  of 
electron  transport  in  biased  GaN/AlGaN  multi-QW  (MQW)  systems,  giving  fi  rst  principal  macro¬ 
scopic  parameters  like  light  and  dark  current  density,  responsivity,  etc.,  (ii)  design,  modelling  and 
experimental  verifi  cation  of  a  GaN/AlGaN  QWIP  for  detection  in  the  near-infrared  spectral  range. 


THEORETICAL  CONSIDERATIONS 

We  consider  P  periods  of  a  MQW  structure  with  an  externally  applied  electric  fi  eld.  The  energy 
spectrum  of  such  a  system  in  principal  is  continuous  which  makes  the  exact  description  of  electron 
behavior  virtually  impossible.  In  order  to  overcome  this  here  we  introduce  two  main  approxima¬ 
tions:  i)  that  the  continuum  can  be  accurately  represented  with  a  number  of  discrete  states  selected 
from  the  continuum  by  an  appropriate  boundary  condition,  ii)  that  the  MQW  structure  is  ideally 
periodic  (and  all  periods  are  equivalent),  which  gives  an  unique  set  of  states  associated  to  a  certain 
period. 

In  such  a  case  the  electron  transport  can  be  simulated  by  implementing  a  rate  equation  ap¬ 
proach,  from  which  the  electron  distribution  over  quantized  continuum  states  can  be  in  steady-state 
calculated  by  solving  a  system  of  linear  equations  as: 
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where  i  +  kN  is  the  /th  state  of  A:th  neighbor  period,  W;  ^  is  the  total  scattering  rate  from  state  i 
into  state  j,  rii  is  the  electron  concentration  of  the  ith  state  and  ^  is  equal  to  1  for  light  and  0 
for  dark  conditions.  The  parameter  C(.)  is  equivalent  to  a  external  perturbation,  in  our  case  the 
concentration  of  electrons  generated  per  second  by  the  applied  optical  fi  eld.  For  known  electron 
distributions  over  continuum  states  one  can  estimate  macroscopic  parameters  of  the  system  like 
current  density,  capture  probability,  responsivity  etc. 

The  current  density  can  be  calculated  by  subtracting  the  current  density  component  due  to 
electrons  scattering  into  the  next  period  from  the  component  due  to  electrons  scattering  back,  as 
usually  used  in  quantum  cascade  laser  simulations  [16].  As  we  are  dealing  with  continuum  states 
this  approach  can  be  intrinsically  instable.  Hence  we  introduce  a  different  method  assuming  that 
average  speed  of  the  electron  scattering  from  the  state  i  to  j  is  given  as 


(vij)  = 


(zj)  -  {Zi) 


(2) 


where  X/,-  is  the  electron  scattering  time  between  states  i  and  j  (X/.-  =  VF  , 2),  z  is  an  average  coordi- 


hj 


nate  of  the  electron  given  as 


{zi)  =  J  \\fi{z)z\\fUz)dz  (3) 

—  CO 

where  \|/,-  is  the  wavefunction  of  the  ith  state  and  z  coordinate  along  the  growth  direction.  Then  the 
current  density  can  be  calculated  as 


p/2-1  N  N  ^  P/2-1  N  N  ^ 


(4) 
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where  Ns  is  total  sheet  electron  density  in  one  period,  L  is  the  length  of  the  period.  This  expression 
deli  nes  both  the  dark  and  the  light  current  dependent  on  the  value  of  the  parameter  When  the 
current  density  is  known,  the  peak  responsivity  can  be  expressed  as 


7(^=1) -7(^  =  0) 
hcj\ 


where  X  is  the  detection  wavelength. 
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A  very  important  issue  is  the  choice  of  wavefunctions  assigned  to  a  certain  period,  as  the 
continuum  states  are  not  confi  ned  within  a  deli  ned  area  and  the  wrong  choice  can  lead  to  the 
existence  of  duplicated  states  and  a  false  estimation  of  the  carrier  dynamics.  In  order  to  overcome 
this  potential  problem  we  deli  ne  states  belonging  to  the  considered  period  as  those  having  the  better 
overlap  integral  with  the  ground  state  of  that  period  than  with  the  ground  states  of  the  neighboring 
periods.  For  this  purpose  we  deli  ne  the  overlap  integral  as: 

oo 

Oij=  j  (6) 

—  CO 

In  order  to  reduce  the  number  of  scattering  rate  processes  necessary  to  calculate  the  electron 
distribution  and  corresponding  current  density  (note  that  the  number  of  total  scattering  rate  pro¬ 
cesses  is  equal  to  N^{2P+  1)  —N),  we  introduce  the  ‘tight-binding’  approximation  assuming  that 
only  a  few  closest  neighbors  interact,  and  set  P  =  2.  Also  we  assumed  that  electron  -longitudinal 
optical  (LO)  phonon  scattering  is  the  main  and  only  scattering  mechanism.  Electron-electron  scat¬ 
tering  was  considered  to  be  negligible  as  the  continuum  states  are  hardly  populated.  The  A:-space 
scattering  rate  averaging  assumed  Fermi-Dirac  distribution  within  each  subband  with  a  unique 
electron  temperature. 

The  energies  and  wavefunctions  were  calculated  by  solving  the  envelope  function  Schrodinger 
equation,  within  the  effective  mass  approximation,  with  the  effects  of  bulk  nonparabolicity  via 
the  energy  dependent  effective  mass  (two  band  Kane  model)  and  many-body  exchange  potential 
included.  The  large  intrinsic  electric  fi  eld  (F  ~  MV/cm)  induced  by  the  nonlinear  piezoelectric 
and  spontaneous  polarizations  is  introduced  by  assuming  periodic  boundary  conditions  over  the 
whole  structure  [17,  18]  as: 

L{Pk-Pj)lk/ek 

ejLk/ek 

k 

where  Pi  =  Pip  -\-  Pp^  is  the  total  polarization  given  as  a  sum  of  piezoelectric  and  spontaneous 
(pyroelectric)  polarization  within  a  layer  of  width  //.  In  calculating  the  alloy  band  gap  we  have 
included  the  effect  of  bowing  [19]. 

The  optical  perturbation  constant,  in  the  steady  state  conditions,  is  modelled  assuming  linear  a 
dependence  on  the  optical  fi  eld  flux  as: 

C(d>,n/,n/)  =A,y$  (8) 

where  Ay  is  fractional  optical  apsorption  on  the  i  f  transition  which  reads: 

/y^ft\  /*cxD 

Ay((o)  =  |My|  L{~h&~}mkj)Fif{kj)dkj  (9) 
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where  COq  is  the  transition  frequency,  ~h  is  the  photon  energy,  c  the  velocity  of  light  in  vacuum,“h  the 
average  refraction  index,  Mif  the  dipole  matrix  element,  kt  the  transversal  (in-plane)  wavevector 
{kx,ky},  Fif  is  the  difference  of  Fermi-Dirac  factors  for  the  two  states,  and  L(.)  the  normalized 
Lorentzian: 


L(  h(D,  h§))kj) 
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where  F  is  the  homogeneous  part  of  the  transition  linewidth  (FWHM),  and  the  transition  energy 
“/j(Q  depends  on  kt,  i.e. 


'h(Q  =  Eof  —  Eoi  + 
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where  Fq/,/  denote  the  state  energies  at  kt  =  0,  and 


(12) 


is  the  difference  of  the  reciprocal  in-plane  (transversal)  electron  effective  masses  my  in  states  i 
and  /.  These  masses  are  unequal  due  to  both  the  nonparabolicity  and  the  position  (material) 
dependence  of  the  band  edge  effective  mass. 

The  model  presented  above  is  based  on  a  fi  rst  principal  calculation  with  only  the  optical  per¬ 
turbation  as  an  input  parameter.  Commonly  used  methods  can  not  deal  with  only  a  few  discrete 
states  in  the  quantum  wells  [13],  which  is  overcome  here  by  the  generality  of  this  approach.  Its 
worth  mentioning  that  a  few  other  authors  applied  quantum  mechanical  scattering  rate  analysis  for 
dark  current  and  capture  probability  calculations  [20,  21],  giving  no  insight  into  electron  transport 
properties. 


NUMERICAL  RESULTS 

In  this  section  we  present  design,  theoretical  modelling  and  experimental  confi  rmation  of  the 
GaN/AlGaN  QWIP.  The  structure  was  designed  for  intersubband  transitions  between  the  lowest 
two  subbands  of  around  250meV  corresponding  to  a  wavelength  of  4/rm  in  the  non-biased  con¬ 
ditions.  Also  the  strain-balancing  requirement  was  enforced  by  using  the  ’zero-strain‘  approach 
[12].  The  MQW  structure  grown  by  MOVPE  [22],  consisted  of  10  QWs  with  60 A  wide  wells 
and  280A  wide  barriers.  The  A1  content  in  the  barrier  layers  was  set  to  40%  giving  the  equiva¬ 
lent  barrier  height  of  0.75eV  (AE^  ~  1 .8eV).  The  Raman  spectroscopy  study  revealed  intersubband 
transitions  at  around  4. 1/rm  [22],  and  confi  rmed  our  choice  of  alloy  parameters  as  well  as  the  appli¬ 
cation  of  our  self-consistent  model  for  the  calculation  of  intersubband  energies  and  corresponding 
wavefunctions. 
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Wavelength  [|am] 

FIG.  1:  Calculated  absorption  profi  les  for  biased  (solid  line)  and  non-biased  (dashed  line)  structure. 

When  an  external  electric  fi  eld  is  applied,  a  shift  of  the  intersubband  transition  energies  occurs 
dependent  on  the  strength  and  polarization  of  the  fi  eld.  Fig.  1  presents  the  calculated  intersubband 
absorption  profi  les  (based  on  the  model  presented  in  the  previous  section)  for  zero  bias  and  an 
applied  reverse  bias  of  20V.  The  estimated  blue  shift  is  around  1/rm  suggesting  a  wide  tunability. 

Fig.  2  shows  the  sample  dark  current  measured  at  77K.  The  TV  curve  is  asymmetric  due  to 


Bias  [V] 

FIG.  2:  Measured  dark  current  vs.  applied  bias  at  77K  for  GaN/AlGaN  QWIP.  Dots  present  theoretical 
predictions. 
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FIG.  3:  Measured  responsivity  vs  wavelength  for  5  (dot-dashed  line),  15  (dashed  line)  and  20V  (solid  line) 
biases. 

the  intrinsic  electric  fi  eld  in  the  layers,  as  the  electrons  can  easier  switch  from  the  quantum  well 
into  continuum,  when  the  externally  applied  electric  fi  elds  is  of  the  same  polarity  as  the  internal 
electric  fi  eld  in  the  barrier.  In  the  working  regime  above  2V  (reverse  bias)  the  value  of  the  dark 
current  is  around  10“^mA.  The  graph  also  presents  theoretical  predictions  for  different  values 
of  applied  bias  and  good  agreement  is  found.  The  slight  discrepancy  is  probably  due  to  other 
leakage  mechanisms  (contacts,  surface  leakage  etc.)  and  are  out  of  the  scope  of  this  work.  We  also 
estimated  the  capture  probability  to  be  around  30%  at  4V  and  99%  at  IV  as  expected.  The  optical 
response  of  the  sample  was  examined  using  the  high  power  source  FELIX  tunable  for  wavelengths 
higher  than  4/rm.  In  Fig.  3  are  presented  responsivity  measurements  as  a  function  of  the  source 
wavelength  for  different  values  of  the  applied  bias  (note  that  the  dips  at  4.2,  6  and  6.5/rm  are  due 
to  atmospheric  absorption).  It  can  be  seen  that  for  20V,  peak  responsivity  occurs  at  around  5 /urn. 
With  the  decrease  of  applied  bias  the  peak  absorption  shifts  towards  shorter  wavelengths  and  as  a 
result  the  optical  response  in  the  range  >  4/rm  decreases.  Fig.  4  shows  the  calculated  responsivity 
as  a  function  of  wavelength  for  experimentally  applied  bias  values  using  the  microscopical  model 
described  in  previous  section.  The  optical  flux  was  assumed  to  be  d>  =10^^cm^^s^^  corresponding 
to  the  average  FELIX  pulse  power.  Again  an  excellent  agreement  with  the  experiment  is  found 
for  both  responsivity  profi  le  and  absolute  values.  The  continuum  states  are  very  weakly  populated 
(lower  than  0.01%  of  the  total  sheet  carrier  density  within  the  period)  even  under  high  power 
excitation  which  opens  the  question  of  an  unexpectedly  high  responsivity,  especially  with  energy 
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Wavelength  [|Ltm] 

FIG.  4:  Calculated  responsivity  vs.  wavelength  for  different  values  of  applied  bias. 

spacings  between  ground  and  continuum  states  above  250  meV.  By  looking  into  scattering  rates  as 
well  as  current  density  components  we  suggest  that  the  high  current  is  due  to  ultra-fast  scattering 
processes  between  pairs  of  continuum  states  separated  by  an  LO  phonon  energy. 


CONCLUSION 

In  conclusion  we  have  presented  a  fi  rst  principals  fully  microscopic  model  for  electron  trans¬ 
port  in  QWIPs  with  an  emphasis  on  structures  with  more  than  one  discrete  state.  The  model  is 
based  on  scattering  rate  equations  with  the  ‘tight-binding’  approximation  applied  by  including  the 
influence  of  two  nearest  neighbors.  Macroscopic  parameters  like  current  density,  responsivity  and 
capture  probability  are  obtained  without  the  need  of  empirical  parameters.  In  order  to  examine  the 
accuracy  of  the  simulation  we  have  applied  it  to  a  GaN/AlGaN  QWIP,  and  found  excellent  agree¬ 
ment  with  the  experiment  for  dark  current  and  responsivity.  The  model  can  generally  be  applied 
to  any  QWIP  device. 
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